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Co-based Heusler compounds are ferromagnetic with a high Curie temperature and a large mag¬ 
netization density, and thus are promising for spintronic applications. In this paper, we perform a 
systematic ab initio study of two principal spin-related phenomena, namely, anomalous Hall effect 
and current spin polarization, in Co 2 Fe-based Heusler compounds Co 2 FeX (X = Al, Ga, In, Si, 
Ge, Sn) in the cubic L2i structure within the density functional theory with the generalized gradi¬ 
ent approximation (GGA). The accurate all-electron full-potential linearized augmented plane-wave 
method is used. First, we find that the spin-polarization of the longitudinal current ( P L ) in Co 2 FeX 
(X = Al, Ga, In, Alo.sSio.s and Sn) is ~100 % even though that of the electronic states at the 
Fermi level ( P D ) is not. Further, the other compounds also have a high current spin polarization 
with P L > 85 %. This indicates that all the Co 2 FeX compounds considered are promising for 
spin-transport devices. Interestingly, P D is negative in Co 2 FeX (X = Si, Ge and Sn), differing in 
sign from the P L as well as that from the transport experiments. Second, the calculated anomalous 
Hall conductivities (AHCs) are moderate, being within 200 S/cm, and agree well with the available 
experiments on highly L2i ordered Co 2 FeSi specimen although they differ significantly from the 
reported experiments on other compounds where the B2 antisite disorders were present. Surpris¬ 
ingly, the AHC in Co 2 FeSi decreases and then changes sign when Si is replaced by Ge and finally by 
Sn. Third, the calculated total magnetic moments agree well with the corresponding experimental 
ones in all the studied compounds except Co 2 FeSi where a difference of 0.3 /rs/f.u. exists. We 
also perform the GGA plus on-site Coulomb interaction U calculations in the GGA+I7 scheme. 
We find that including the U affects the calculated total magnetic moment, spin polarization and 
AHC significantly, and in most cases, unfortunately, results in a disagreement with the available 
experimental results. All these interesting findings are discussed in terms of the underlying band 
structures. 

PACS numbers: 71.20.Be, 72.25.Ba, 75.47.Np, 75.70.Tj 


I. INTRODUCTION 

Most Co-based Heusler compounds in the cubic L2i 
structure are ferromagnetic with a high Curie tempera¬ 
ture and a large saturation magnetization. 1 Furthermore, 
many of them were predicted to be half-metallic 2-5 and 
hence are of particular interest for spintronics. Therefore, 
the electronic band structure and magnetic properties of 
the Co-based Heusler compounds have been intensively 
investigated both theoretically and experimentally in re¬ 
cent years 1-5 . For example, the total magnetic moments 
of these materials were found to follow the Slater-Pauling 
type behavior and the mechanism was explained in terms 
of the calculated electronic structures 2 . The Curie tem¬ 
peratures of Co-based Heulser compounds were also de¬ 
termined from ab initio theoretical calculations and the 
trends were related to the electronic structures 3 . 

Half-metallic ferromagnets are characterized by the co¬ 
existence of metallic behavior for one spin channel and 
insulating behavior for the other, and their electronic 
density of states at the Fermi level is completely spin 
polarized. Thus, they could in principle offer a fully 
spin-polarized current and are useful for spin electronic 
devices. The possible half-metallicity of the Co-based 
Heusler compounds has been intensively investigated 


experimentally 6-15 and theoretically 2-5,7,10 . In partic¬ 
ular, many point-contact Andreev reflection (PCAR) 
experiments 6-10 have been carried out on Co 2 FeSi and its 
current spin polarization ( P L ) was found to vary from 45 
% to 60 %, depending on the substrate and the quality 
of the contact. A higher P L of ~80 % was reported in 
a nonlocal spin valve (NLSV) experiment 11 on Co 2 FeSi. 
The spin Hall effect experiment 8 showed that the P L 
of Co 2 FeSi is positive. However, the measured positive 
current spin polarization is at odds with the predictions 
of the negative spin polarization of the electronic states 
at the Fermi level (static spin polarization) ( P D ) from 
the ab initio calculations 7,16 with the local density ap¬ 
proximation (LDA) or generalized gradient approxima¬ 
tion (GGA) 17 . Furthermore, the calculated spin mag¬ 
netic moment (m| ot ) was found to differ by nearly 10 % 
from the measured total magnetic moment ( m toi ) of ^6 
PLb/ f.u. 16 It was argued that the Co-based Heusler com¬ 
pounds are strongly correlated systems and hence should 
be better described by, e.g., the LDA/GGA plus onsite 
Coulomb Repulsion ( U ) (LDA/GGA+U) approach 18 . 
The GGA+t/ calculations 16 indeed would give rise to a 
?n* ot of 6 ^b/ f-u. and a positive P D of 100 % (i.e., be¬ 
ing half-metallic). However, the calculated P D is much 
larger than the measured spin polarization. Nevertheless, 
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as pointed out recently in Ref. 7, the spin polarization 
measured in transport experiments such as PCAR and 
NLSV experiments should be compared with the theo¬ 
retical current spin polarization ( P L ) rather than static 
spin polarization ( P D ). Therefore, one of the principal 
purposes of this work is to understand the measured spin 
polarization in all the Co 2 FeX compounds by performing 
a systematic ab initio GGA study of both the current and 
static spin polarizations as well as the total magnetic 
moment in these compounds. Indeed, we find that the 
calculated P L values for the Co 2 FeX compounds agree 
with the available experimental results not only in sign 
but also in magnitude, while the calculated P D values 
for Co 2 FeX (X = Si, Ge and Sn) are wrong even in sign. 

The anomalous Hall effect (AHE), discovered in 
1881 by Hall 19 , is an archetypal spin-related transport 
phenomenon and hence has recently received renewed 
attention 20 . Indeed, many ab initio studies on the 
AHE in elemental ferromagnets 21 24 and intermetallic 
alloys 25,26 have recently been reported. However, ab ini¬ 
tio investigations into the AHE in the Heusler compounds 
have been few 5,27,28 . Interestingly, Co 2 MnX (X = Al, Ga 
and In) were recently predicted 5 to have a large intrin¬ 
sic anomalous Hall conductivity (AHC) in the order of 
~1000 S/cm, and thus could find applications in magne¬ 
tization sensors 29 . Therefore, another principal purpose 
of this work is to understand the AHE in all the Co 2 Fe- 
based Heusler compounds and the results may help ex¬ 
perimental search for the Heusler compounds with large 
AHE for applications. Furthermore, by comparison of the 
calculated AHC as well as the current spin polarization 
and total magnetic moment with the measured ones, one 
could have a comprehensive assessment of whether or not 
the Co 2 FeX compounds are strongly correlated systems 
that would require the GGA+I7 approach. 

In nonmagnetic materials where the numbers of the 
spin-up and spin-down electrons are equal, the opposite 
transverse currents caused by the applied electric field 
would result in a pure spin current, and this is known 
as the intrinsic spin Hall effect (SHE ). 30 The pure spin 
current is dissipationless 30 and is thus important for the 
development of low energy-consumption nanoscale spin- 
tronic devices 31 . We note that high spin-polarization 
( P L ) of the charge current (Ic) from the electrode is es¬ 
sential for large giant magnetoresistance (GMR ) 32,33 and 
tunneling magnetoresistance (TMR) 34,35 . However, since 
the current-induced magnetization switching results from 
the transfer of spin angular momentum from the current 
carriers to the magnet 36 , large spin current (Is) would 
be needed for the operation of the spin-torque switching- 
based nanodevices 36,37 , i.e., a large ratio of spin current 
to charge current [77 = \(2e/H)Is / Ic\, would be crucial. 
For ordinary charge currents, this ratio 77 varies from 
0.0 (spin unpolarized current) to 1.0 (fully spin polar¬ 
ized current). Interestingly, 77 can be larger than 1.0 for 
the Hall currents and is 00 for pure spin current. Fasci¬ 
natingly, spin-torque switching of ferromagnets driven by 
pure spin current from large SHE in tantalum has been 


recently reported 31 . Therefore, it might be advantageous 
to use the Hall current from ferromagnets for magneto- 
electronic devices, rather than the longitudinal current. 
Another purpose of this work is therefore to investigate 
the nature and spin-polarization of the Hall current in the 
Co 2 Fe-based Heusler compounds for possible spintronic 
applications. 

The rest of this paper is organized as follows. In the 
next section, we briefly describe the theories of the intrin¬ 
sic anomalous and spin Hall conductivities as well as Hall 
and longitudinal current spin polarizations. We will also 
introduce the full-potential relativistic band theoretical 
method used and give the computional details. In Sec. 
Ill, the calculated magnetic moments, intrinsic Hall con¬ 
ductivities and current spin polarizations of all the stud¬ 
ied Co 2 FeX compounds will be reported in subsections 
III (A), (B) and (C), respectively. The theoretical results 
will be analyzed in terms of the underlying band struc¬ 
tures and also compared with the available experimental 
ones. In subsection III (D), the results from the GGA+U 
calculations will be presented to examine the effect of in¬ 
cluding the semi-empirical on-site Coulomb interaction 
on the calculated physical properties of the Co 2 FeX com¬ 
pounds. Finally, conclusions drawn from this work will 
be given in Sec. IV. 


II. THEORY AND COMPUTATIONAL 
METHOD 

We first perform the self-consistent electronic structure 
calculations for the Co 2 FeX compounds within the den¬ 
sity functional theory with the GGA for the exchange cor¬ 
relation potential 17 . Since all the intrinsic Hall effects are 
caused by the relativistic electron spin-orbit interaction, 
the spin-orbit coupling (SOC) is included in the present 
ab initio calculations. We use the highly accurate full- 
potential linearized augmented plane wave (FLAPW) 
method, as implemented in the WIEN2K code 38 . The 
wave function, charge density, and potential were ex¬ 
panded in terms of the spherical harmonics inside the 
muffin-tin spheres and the cutoff angular moment ( L max ) 
used is 10, 6 and 6 , respectively. The wave function out¬ 
side the mufRn-tin spheres is expanded in terms of the 
augmented plane waves (APWs) and a large number of 
APWs (about 70 APWs per atom, i. e., the maximum 
size of the crystal momentum K max = 8/R m t) were in¬ 
cluded in the present calculations. The improved tetrahe¬ 
dron method is used for the Brillouin-zone integration 39 . 
To obtain accurate ground state charge density as well 
as spin and orbital magnetic moments, a fine 27x27x27 
grid with 1470 fc-points in the irreducible Brillouin zone 
wedge (IBZW) is used. The self-consistent cycles were 
terminated when the integrated charge density variation 
became less than 10~ 5 e. 

We consider the Co 2 FeX Heusler compounds in the 
fully ordered cubic L2i structure. The available exper¬ 
imental lattice constants 40 are used for all the consid- 
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ered Co 2 FeX (X = Al, Ga, In, Si, Ge, Sn) Heusler alloys 
except Co 2 FeAl 0 . 5 Sio. 5 , Co 2 FeIn and Co 2 FeSn, as listed 
in Table I. Since the experimental lattice constant for 
Co 2 FeIn is not available, we determine the lattice con¬ 
stant for Co 2 FeIn theoretically, also by using the FLAPW 
method, as described in the preceding paragraph. We 
also study the L2i Co 2 FeAlo. 5 Sio .5 alloy and model it by 
the virtual crystal approximation (VGA), i.e., the Al/Si 
site is occupied by a virtual atom with the atom number 
Z = 0 . 5 . 214 ; + 0.52Tgj, where Zai and Zgi are the Al and 
Si atomic numbers, respectively. The lattice constant of 
5.689 A, which is the average of the experimental lattice 
constants of Co 2 FeAl (5.737 A) and Co 2 FeSi (5.640 A), 
is used for Co 2 FeAlo. 5 Sio .5 because the lattice constant of 
the Co 2 FeAl x Sii_ x alloy was reported to depend linearly 
on the Al concentration (a:) 45 . Note that in fact we have 
determined the lattice constants theoretically also for the 
Co 2 FeX (X = Al, Ga, Si, Ge) compounds. The theoreti¬ 
cal lattice constants for these compounds differ from the 
experimental values by less than 1 %. As a result, the 
physical properties of these compounds calculated using 
the experimental and theoretical lattice constants differ 
only slightly. Therefore, for simplicity, we present only 
the physical properties of these compounds calculated 
using the experimental lattice constants in the next sec¬ 
tion. However, the theoretical lattice constant (6.013 A) 
of Co 2 FeSn is 2.4 % larger than the experimental one 
(5.87 A ) 46 perhaps because the prepared Co 2 FeSn films 
contained only a low degree of L2i order. Therefore, the 
theoretical lattice constant is used for Co 2 FeSn (Table I). 


A. Anomalous and spin Hall conductivities 

The intrinsic anomalous and spin Hall conductivi¬ 
ties of a solid can be evaluated by using the Kubo 
formalism. 21,47,48 Here we first calculate the imaginary 
part of the off-diagonal elements of the optical conduc¬ 
tivity. Then we obtain the real part of the off-diagonal 
elements from the imaginary part by a Kramers-Kronig 
transformation. The intrinsic AHC (crA) is the static 
limit of the off-diagonal element of the optical conductiv¬ 
ity ctxy (w = 0). If we now replace the charge current op¬ 
erator —ev with the spin current operator (fi./4){X z ,v} 
and repeat the calculation 4 ', we will obtain the intrinsic 
spin Hall conductivity (SHC) {cr xy ). We note in passing 
that alternatively, one could also calculate a xy {a xy ) by 
an integration of the (spin) Berry curvature over the Bril- 
louin zone 21 ’ 49,50 . Nevertheless, the two methods were 
found to be numerically equivalent 21,49,50 . 

A dense k point mesh would be needed for obtaining 
accurate AHC and SHC 21,48 . Therefore, we use several 
fine /c-point meshes with the finest /e-point mesh being 
58x58x58 which has 8125 /c-points in the IBZW. We 
calculate the AHC and SHC as a function of the number 
( Nk ) of fc-points in the first Brillouin zone. The calcu¬ 
lated AHC (&xy) and SHC (crf y ) versus the inverse of the 
Nk are then plotted and fitted to a polynomial to get the 


converged theoretical a xy and crf y (i.e., the extrapolated 
value at Nk = 00 ) (see Refs. 23 and 24). Furthermore, 
to ensure that the (w = 0) via the Kramers-Kronig 
transformation is accurate, the energy bands up to 5.5 
R.y are included in the calculation of er X j/(w). 


B. Current spin polarization 


The spin polarization of a magnetic material is usually 
described in terms of the spin-decomposed densities of 
states (DOSs) at the Fermi level ( Ep) as follows 


P D = 


N t (E F ) - N^Ep) 
N t (E F ) + Ni(E f ) ’ 


( 1 ) 


where N^(E F ) and N^(E F ) are the spin-up and spin- 
down DOSs at the Ep, respectively. This static spin 
polarization P D would then vary from -1.0 to 1.0 only. 
For the half-metallic materials, P D equals to either -1.0 
or 1.0. As mentioned above, the spin polarization P D de¬ 
fined by Eq. (1) is not necessarily the spin polarization of 
the transport currents measured in experiments. Indeed, 
the spin-polarizations measured by using different exper¬ 
imental techniques could differ significantly 51 " 54 . From 
the viewpoint of spintronic applications, only the current 
spin polarization instead of the P D , counts. 

Therefore, in this work, we further calculate the spin 
polarization of both the longitudinal and Hall currents, 
as described below. Here, we calculate the longitudinal 
electric conductivities (oq,oq) for spin-up and spin-down 
electrons divided by the corresponding Drude relaxation 
times (7q,rq) (i.e., oq/rq, oq/rq) within the semi-classical 
Boltzmann transport theory, as implemented in Boltz- 
Trap code 55 . In the present calculations, the relaxation 
time is assumed to be independent of energy, fc-point and 
spin direction (i.e., rq = tj. = r). Consequently, we can 
obtain the longitudinal current spin polarization P L from 


p L tj 7 (7^ — 

(jt T cjI crl/r + cj 4/r 


( 2 ) 


The underlying scalar-relativistic band structures are cal¬ 
culated by using a fine 36x36x36 mesh with 3349 k- 
points in the IBZW. 

The spin polarization P H of the Hall current may be 
written as 5,24 


cj h t _ a H q 
pH _ xy w xy 


H t 


Hi 


(3) 


where (T^j and <J x J r are the spin-up and spin-down Hall 
conductivities, respectively. The a and can be 
obtained from the calculated AHC and SHC via the 
relations 50 


A _ H t , Hi 
(T — (J (T 

xy xy 1 xy 


(4) 


_9 - a s =a m -a Hi 
Ji xy x y xy ’ 


( 5 ) 
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Note that, the absolute value of P H can be greater than 
1.0 because the spin-decomposed Hall currents can go 
either right (positive) or left (negative). In the non¬ 
magnetic materials, the charge Hall current is zero, and 
hence, = —cr^y results in P H = oo. Clearly, in the 
case of Hall currents, the ratio of the spin current to 
charge current 77 = \{2e/h)a^ y /a^ y \ = \P H |. 

III. RESULTS AND DISCUSSION 
A. Magnetic moments and band structure 

Let us first examine the calculated magnetic properties 
and band structures near the Fermi level of the consid¬ 
ered Co 2 FeX Heusler alloys. Since the electronic struc¬ 
ture and magnetism in the full Heusler compounds have 
been extensively studied (see, e.g., Refs. 2-4 and refer¬ 
ences therein), here we focus on only the salient features 
which may be related to the anomalous and spin Hall ef¬ 
fects as well as current spin polarizations to be presented 
in the next subsections. The calculated total magnetic 
moment, total spin magnetic moment, local spin and 
orbital magnetic moments as well as spin-decomposed 
DOSs at Ep of all the considered Co 2 FeX Heusler alloys 
are listed in Table I, together with the available experi¬ 
mental total magnetic moments for comparison. The to¬ 
tal and site decomposed DOSs of three selected Heusler 
compounds Co 2 FeAl, Co 2 FeAlo. 5 Sio .5 and Co 2 FeSi are 
displayed in Fig. 1. The scalar relativistic band struc¬ 
tures of Co 2 FeAl and Co 2 FeSi are shown in Figs. 2(a) 
and 2(c), respectively. 

Interestingly, the studied Heusler alloys could be sep¬ 
arated into two groups according to the calculated DOS 
at Ep. In one group, including Co 2 FeAl, Co 2 FeGa, 
Co 2 FeAl 0 . 5 Sio .5 and Co 2 FeIn, the majority spin state 
dominates. In the other group, including Co 2 FeSi, 
Co 2 FeGe and Co 2 FeSn, the minority spin state domi¬ 
nates (see Table I). Therefore, the calculated spin polar¬ 
ization ( P D ) for the first group is positive while that of 
the second group is negative (see Table II). It is clear 
from Fig. 1(a) that there is a band gap near the Ep 
for the minority spin channel in Co 2 FeAl and adding one 
valance electron can be approximatively treated as rais¬ 
ing the Ep by ~0.4 eV. This Ep shift is nearly equal 
to that from Co 2 MnAl to Co 2 MnSi in Ref. 5 where the 
calculated P D for both compounds, however, is positive. 
This is because the minority gap here is small, being ~ 
0.2 eV, while that in C^MnAl is much larger, being ~ 
0.8 eV. 

Figure 1 and Table I show that from the view point of 
the calculated band structures, all the considered Heusler 
compounds are not half-metallic, although Co 2 FeAl is 
nearly a half-metal because its Ep just touches the 
bottom of the minority spin conduction band [Fig. 
2(a)]. Previous GGA calculations 58 also predicted that 
of Co 2 FeAl, Co 2 FeAlo. 5 Sio .5 and Co 2 FeSi are not half- 
metallic. The DOS spectra for Co 2 FeAl, C 02 FeAlo. 5 Sio. 5 j 
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FIG. 1. (color online) Total and site decomposed density of 
states (DOS) for (a) Co 2 FeAl, (b) Co2FeAlo.5Sio-5 and (c) 
Co 2 FeSi. The Fermi level is at zero. 


and Co 2 FeSi are similar and differ only in the location of 
Ep (see Fig. 1). The DOS spectra of Co 2 FeGa (Co 2 FeIn) 
and Co 2 FeGe (Co 2 FeSn) (not shown here) also look sim¬ 
ilar, except the location of Ep. 

Table I indicates that among the considered Heusler 
compounds, only the total spin magnetic moment m* ot 
of Co 2 FeAl and Co 2 FeGa almost satisfies the so-called 
generalized Slater-Pauling rule ?n‘ ot = n v —24 where n v is 
the number of valence electrons 2 . This may be expected 
because none of these compounds is predicted to be half- 
metallic here and only Co 2 FeAl is nearly half-metallic. 
Table I also suggests that the calculated m tot agrees well 
with the available measured one for all the considered 
compounds except Co 2 FeSi. The discrepancy between 
the calculated and experimental m tot is ~0.3 /rs/f.u. for 
Co 2 FeSi but is about 0.1 ^s/f.u. or less for all the other 
compounds (Table I). 


B. Anomalous and spin Hall conductivities 

The calculated anomalous Hall conductivity cr£ y and 
spin Hall conductivity a^ y for all the studied compounds 
are listed in Table II. We notice that compared with the 
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TABLE I. Calculated total spin magnetic moment (m“) (/is/f.u.), atomic spin (m s ) and orbital ( m 0 ) magnetic moments 
(/is/atom) as well as spin-decomposed density of states at the Fermi level [N^(Ep), N^(Ef)] (states/eV/f.u.) of all the 
considered Co 2 FeX Heusler compounds together with the lattice constants a (A) used. The available experimental mag¬ 
netic moments 41-44 (Exp.) are also listed for comparison with the calculated total magnetic moments ( m tot ) (/ts/f.u.). In 
Co 2 FeAlo. 5 Sio. 5 , listed in the bracket is the spin magnetic moment of Si. The orbital magnetic moments for the non-transition 
metal atoms (m*) are less than 0.0001 /is/atom and hence not listed here. The theoretical total magnetic moment ( m tot ) is 
given by m tot = + 2which should be compared with the experimental magnetic moment. 


Co 2 FeX 

a 

m tot 

ml oi 

m c :° 

mf 

mf 

m^° 

m o B 

N^{E F ) 

N^(E f ) 

Co 2 FeAl 

5.737 a GGA 
Exp. 

5.123 

4.96 

4.993 

1.229 

2.788 

-0.064 

0.041 

0.048 

0.862 

0.059 

Co 2 FeGa 

5.751 a GGA 
Exp. 

5.149 

5.13 

5.016 

1.206 

2.811 

-0.047 

0.041 

0.051 

0.885 

0.189 

CoaFeln 

5.990 GGA 

5.308 

5.143 

1.250 

2.885 

-0.046 

0.052 

0.061 

0.859 

0.575 

C 02 Fe Alo. 5 Sio.5 

5.689 GGA 

5.523 

5.376 

1.338 

2.683 

-0.037 

0.052 

0.043 

0.755 

0.399 

CoaFeSi 

5.640 a GGA 
Exp. 

5.688 

5.97 

5.541 

1.388 

2.848 

-0.002 

0.040 

0.067 

0.714 

2.476 

CoaFeGe 

5.743 a GGA 
Exp. 

5.854 
5.90, 5.74 

5.693 

1.422 

2.917 

0.012 

0.046 

0.069 

0.785 

2.288 

CoaFeSn 

6.013 GGA 

5.994 

5.797 

1.445 

3.021 

-0.005 

0.060 

0.079 

0.712 

2.457 


a Experimental lattice constants 40 . 


TABLE II. Calculated anomalous [a^y (S/cm)] and spin [a^ y (fiS/e cm)] Hall conductivities, spin-decomposed Hall conductivities 
{<Txy, o'xy') (S/cm), Hall ( P H ) and longitudinal ( P L ) current spin polarizations (%) as well as spin polarization of the electronic 
states at the Fermi level P D (%) of all the considered Heusler compounds Co 2 FeX. The available experimental spin polarization 
and scattering-independent part (6) 56 of the a£ y are also listed for comparison. Note that b contains both the intrinsic 
contribution (a^y) calculated here and also the extrinsic side-jump contribution (<t/ 1 j) _ s - 7 ). 


Co 2 FeX 


A 1 

&xy 

s 

v xy 

a* t 

xy 

a Hl 

P H 

pL 

pu 

CoaFeAl 

GGA 

39 

35 

-16 

55 

-180 

100 

87 


Exp. 

320~360 a 





56 b 


CoaFeGa 

GGA 

181 

56 

35 

147 

-62 

98 

65 


Exp. 






57 c 


CoaFeln 

GGA 

102 

56 

-5 

107 

-110 

92 

20 

Goa Fe Alo. 5 Sio. 5 

GGA 

124 

74 

-12 

136 

-119 

92 

31 


Exp. 

-100~50 s , 352 h 





60 e 


CoaFeSi 

GGA 

189 

24 

71 

119 

-25 

86 

-55 


Exp. 

163 d ,300~600 a 





45~60 e , 80 f 


CoaFeGe 

GGA 

119 

-29 

89 

31 

49 

89 

-49 


Exp. 






59 c 


CoaFeSn 

GGA 

-78 

-24 

-15 

-63 

-62 

93 

-55 


a Experimental b values from sputtered films with the B2 structure 13 . 
b Point-contact Andreev reflection experiments 12 . 

c Point-contact Andreev reflection experiments on C 02 FeGa„ ; Gei - x in the L2i/B2 mixed structure 15 . 
d Experimental b value from Co 2 FeSi single crystals with the L2i structure®. 
e Point-contact Andreev reflection experiments 7 ’ 0,10,45 . 

1 Nonlocal spin-valve experiment 11 . 

g Experimental b values from sputtered Co 2 FeAlo. 4 Sio .6 films with the B2 structure 57 . 

1 Experimental b value from sputterd ultrathin Co 2 FeAlo. 5 Sio .5 film with the B2 structure 14 . 


AHC of Fe metal 21 and also C 02 M 11 X (X = Al, Ga and 
In) 5 , the u^ y of the present Heusler compounds is mod¬ 
erate in magnitude, being within 200 S/cm (Table II). 
In fact, the AHC of Co 2 FeX (X = Al, Ga and In) (Ta¬ 
ble II) is about one order of magnitude smaller than the 
corresponding Co 2 MnX (X = Al, Ga and In) (see Table 


II in Ref. 5). This can be explained in terms of the cal¬ 
culated band structure and also cr/L as a function of E F 
in Co 2 FeAl (Figs. 2a and 2b). Figures 2a and 2c show 
that the spin-up bands near Ep are the highly disper¬ 
sive Co spd, Fe spd and Al sp hybridized bands while 
Ep nearly falls within the spin-down band gap. Conse- 
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quently, the <j£ y i s rather small (being ~35 S/cm) (see 
Fig. 2b). However, when Ep is lowered to below -0.8 eV, 
a^ y increases dramatically to the values of ~1000 S/cm 
(Fig. 2b). These large a^ y values come mainly from the 
spin-up Co d dominant bands in this energy range (Figs. 
2a and 2c). In C 02 M 11 AI, the corresponding spin-up Co 
d dominant bands are higher in energy, and the Ef is 
lower because Co 2 MnAl has one fewer valence electron 
than Co 2 FeAl. As a result, the Ep sits on the Co d dom¬ 
inant (J^y peak in Co 2 MnAl and thus Co 2 MnAl has a 
much larger <j£ (being ~ 1300 S/cm ). 5 This interesting 
finding suggests a way to chemical composition tuning of 
the AHC in Co 2 Mn 1 _ a ,Fe a: X (X = Al, Ga and In) alloys. 

Interestingly, for the Co 2 FeX (X = Si, Ge and Sn) 
compounds, the AHC gets reduced when Si is replaced 
by Ge and changes sign when Ge is further substituted 
by Sn. Nevertheless, the calculated band structures for 
the Co 2 FeX (X = Si, Ge and Sn) compounds look very 
similar, especilly in the vinicity of Ep. Thus, there is no 
obvious explanation for this interesting evolution. Ta¬ 
ble II indicates that the of Co 2 FeSi is about five 
times larger than that of Co 2 FeAl. This could be at¬ 
tributed to the band filling effect. Figure 2 shows that in 
Co 2 FeSi, due to the additional one valence electron, the 
Ep is raised to the bottom of the Co/Fe d(e g ) dominant 
bands where a^ y is large (Fig. 2b), thus resulting in a 
much larger a^y 

Several AHE experiments on the Co 2 FeX compounds 
and their alloys have been carried out . 6,13,14,57 The de¬ 
rived AHC values ( 6) 56 for C 02 FCAI, Co 2 FeAlo. 5 Sio .5 and 
Co 2 FeSi are listed in Table II. However, quantitative 
comparison of the present theoretical calculations with 
the experimental results is difficult, because all the sam¬ 
ples used in the experiments except Co 2 FeSi are in the 
B2 structure with antisite disorders (Table II). The de¬ 
duced AHC values depend strongly on the substrates 
used and annealing temperatures which control the de¬ 
gree of the B 2 antisite disorders and also the defect 
concentrations . 13,57 Nevertheless, Table II shows that the 
calculated cr^ y of Co 2 FeSi is in good agreement with the 
experimental result from the single crystal sample 6 , indi¬ 
cating the intrinsic AHC a^ y dominates in Co 2 FeSi single 
crystals with the L2i structure. In contrast, the theoreti¬ 
cal a^y of Co 2 FeAl is one order of magnitude smaller than 
the b derived from the experiment 13 . For Co 2 FeAl 0 . 5 Sio. 5 , 
the b values from two different experiments 14,57 are very 
different, suggesting the important influences of the B2 
antisite disorder and also the substrate. We could at¬ 
tribute the pronounced discrepancies between the theo¬ 
retical (intrinsic) (a^y) and experimental b values 13 of 
C 02 FCAI and Co 2 FeSi to a significant contribution from 
the impurity side-jump scattering as well as the struc¬ 
tural difference. However, recent ab initio calculations 2 ' 
for Co 2 CrAl and Co 2 MnAl indicated that the B2 antisite 
disorders tend to significantly reduce the intrinsic AHC 
a^y. Therefore, in the experiments 13 on Co 2 FeAl and 
Co 2 FeSi, side-jump mechanism could dominate and thus 


result in a much larger b than a£. 

Table II indicates that the calculated a (L in Co 2 FeIn 

•by 

and Co 2 FeAlo. 5 Sio .5 is about half of the and their 
Hall current spin polarization ( P H ) is nearly 100 %. In a 
half-metal, the charge current would flow only in one spin 
channel and no charge current in the other spin channel, 
thus resulting in a^ y being twice as large as a^ y . There¬ 
fore, Co 2 FeIn and Co 2 FeAl 0 . 5 Sio .5 may be called anoma¬ 
lous Hall half-metals 5 , even though their electronic states 
near Ep are far from fully spin-polarized (see P D in Ta¬ 
ble II). Finally, we note that the ratio of spin current to 
charge current for the Hall current (77 = \P H |) in Co 2 FeAl 
is large with 77 > 150 %. 

C. Current spin polarizations 

The calculated spin polarizations of Hall ( P H ) and lon¬ 
gitudinal ( P L ) currents as well as electronic states at Ep 
( P D ) for all the Heusler componds considered here are 
listed in Table II. Also listed in Table II are the spin- 
decomposed Hall conductivities and (J^y ) obtained 
using Eqs. (4) and (5). Remarkably, Table II shows that 
the calculated P L is nearly 100 % in Co 2 FeAl, Co 2 FeGa 
and Co 2 FeAl 0 . 5 Sio .5 even though their P D is significantly 
smaller than 100 %. This finding, therefore, indicates 
that these Heusler compounds are half-metallic from the 
viewpoint of charge transport, even though their elec¬ 
tronic band structures are not. All the other compounds 
also have a high current spin polarization with P L > 85 
%. Therefore, all the Heusler compounds considered here 
may find valuable applications in spintronic devices. In¬ 
terestingly, Table II also demonstrates that the P L and 
P D in Co 2 FeSi, Co 2 FeGe and Co 2 FeSn could even have 
opposite signs. The calculated current spin polarization 
P L in Co 2 FeSi and Co 2 FeGe is positive, being in good 
agreement with recent spin Hall effect experiments 8 . In 
contrast, the static spin polarization ( P D ) differs from 
the experimental spin polarization even in sign (Table 
II). This clearly urges one to compare the measured spin 
polarization from transport experiments to the theoreti¬ 
cal current spin polarization rather than the static spin 
polarization which has often been done in the past. 

The interesting finding that the P L and P D in 
Co 2 FeSi, Co 2 FeGe and Co 2 FeSn differ in sign, could be 
explained in terms of the calculated band structures. Fig¬ 
ure 2 (c) indicates that in Co 2 FeSi, for the spin-up chan¬ 
nel, the Ep cuts through the highly dispersive Co/Fe 
spd and Si sp hybridized bands. On the other hand, 
for the spin-down channel, the Ep is located at the bot¬ 
tom of the Co/Fe d(e g ) dominated bands. Consequently, 
the spin-down DOS at Ep is higher than the spin-up 
DOS (see Fig. lc and Table I), giving rise to the neg¬ 
ative value of P D . From transport viewpoint, however, 
the spin-down Co/Fe d{e g ) dominated bands which are 
narrow (Fig. 2c), would have large effective masses and 
small Fermi velocities, thereby contributing little to the 
charge current. On the other hand, the spin-up Co/Fe 
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FIG. 2. (color online) Scalar relativistic band structure [(a) 
and (c)], and anomalous Hall conductivity (a£ y ) [(b) and (d)[ 
for CogFeAl, and Co 2 FeSi, respectively. The Fermi energy is 
at zero. 


spd and Si sp hybridized bands which are highly disper¬ 
sive, would have large Fermi velocities and small effective 
masses, thus providing the dominant contribution to the 
charge current. Therefore, the current spin polarization 
is positive, being in good agreement with the experiments 
(Table II). 

Many experiments 7,9,10 ’ 12 ’ 15 ’ 45 especially PCAR mea¬ 
surements on Co 2 Fe-based Heusler alloys have been car¬ 
ried out to determine their spin polarization which is a 
key factor for their spintronic applications. Majority of 
these experiments were focused on Co 2 FeSi mainly be¬ 
cause highly L2i ordered Co 2 FeSi samples could be fab¬ 
ricated. However, the P L values derived from PCAR 
experiments on Co 2 FeSi vary significantly from 45 % to 
60% (Table II), depending on the quality of the sam¬ 
ples. This could be expected because the spin polar¬ 
ization determined by a PCAR experiment depends not 
only on the degree of the ordering and the defects in the 
sample but also on the quality of the contact and the 
substrate 59 . Nevertheless, the theoretical P L value of 86 
% agrees rather well with the experimental value of 80% 
from the nonlocal spin-valve experiment 11 on highly L2i- 
ordered specimens. However, the measured P L values for 
Co 2 FeAl, Co 2 FeGa and Co 2 FeAl 0 . 5 Sio .5 are around 60 %, 
which is far from the predicted P L value of ~100 % for 
these compounds (Table II). These significant discrepan¬ 
cies may reflect the fact that the samples used in the 


experiments 10,12,15 had a high degree of the B2 antisite 
disorders. 


D. Effects of on-site Coulomb interaction 

To examine the effect of on-site Coulomb interac¬ 
tion, we further perform the calculations in the GGA+C/ 
scheme 18 . The on-site Coulomb repulsion [/(exchange 
interaction J) used are 2.82 (0.9) and 2.6 (0.8) eV for Co 
and Fe, respectively, which are widely used for Co 2 -based 
Heusler compounds. 58 The results from these GGA+C/ 
calculations are compared with those of the GGA calcu¬ 
lations in Table III. We notice that the total spin mag¬ 
netic moments (m*°*) from the GGA and GGA+C/ calcu¬ 
lations are almost identical in all the Heusler compounds 
except Co 2 FeSi and Co 2 FeGe. This may be expected 
since the GGA m* ot is already nearly saturated in these 
compounds. Including the on-site Coulomb interaction 
increases the ?n* ot in Co 2 FeSi and Co 2 FeGe to the satu¬ 
ration values. Note that the measured magnetic moments 
should be compared with the calculated total magnetic 
moments ( m tot ) instead of total spin magnetic moments 
(ro* ot ) in Tables I and III. The m tot contains both the 
to*°* and the total orbital magnetic moment which can¬ 
not be neglected in the Heusler compounds studied here 
because the orbital magnetic moments on the Fe and Co 
atoms are rather significant (Table I). Tables I and III 
together show that including the on-site Coulomb inter¬ 
action actually increases the small discrepancies between 
the experiments and the GGA calculations found in all 
the Heusler compounds except Co 2 FeSi where the dif¬ 
ference of 0.3 ps/ f-u. is reduced slightly to 0.2 ps per 
formula unit. 

Table III shows that the on-site Coulomb interaction 
has a pronounced effect on the spin polarization. First, 
the current spin polarization ( P L ) for all the studied com¬ 
pounds is now 100 % from the GGA+C/ calculations. 
Second, the static spin polarization ( P D ) approaches to 
100% for all the compounds except Co 2 FeGa. Therefore, 
these Heusler compounds become half-metals in terms 
of both the band structure and current spin polariza¬ 
tion. This may be expected because the main effect of 
the on-site Coulomb interaction is to raise the spin-down 
Fe and Co d-dominant conduction bands. Consequently, 
if sufficiently large U values are used, the spin-down Fe 
and Co d-dominant conduction bands will move to above 
E p. This will open a gap in the spin-down channel 
and thus give rise to zero spin-down DOS at Ep (Ta¬ 
ble III). However, the spin-down GGA+C/ band gap is as 
large as 0.9 eV in Co 2 FeSi, for example, being nearly 100 
times larger than the measured one 6 . Interestingly, in¬ 
cluding Coulomb U changes the spin polarization P D in 
Co 2 FeSi and Co 2 FeGe from negative to positive (Table 
III). However, it should be emphasized that the mech¬ 
anism of the spin polarization sign change here is very 
different from the sign difference between the P D and 
P L in the GGA calculations. Nevertheless, whether the 









































TABLE III. Total magnetic moment ( m tot ), total spin magnetic moment (ml 01 ) (/is/f-U.), spin-decomposed density of states 
at the Fermi level | N^(E F ), N^(E f )] (states/eV/f.u.), spin polarization of the electronic states at the Fermi level P D (%), 
longitudinal current polarization P L (%, anomalous [a^ y (S/cm)] and spin [a^ y ( hS/e cm)] Hall conductivities and Hall current 
spin polarization P H (%) of the Co 2 FeX Heusler compounds from both the GGA and GGA+U calculations. The on-site 
Coulomb (exchange) interaction U (J) for Co and Fe used are 2.82 (0.9) eV and 2.6 (0.8) eV, respectively. 


Co 2 FeX 

m tot 

ml ot 

N r (E F ) 

N±(E f ) 

pu 

pL 

V xy 

®xy 

pH 

Co 2 FeAl 

GGA 5.123 

4.993 

0.862 

0.059 

87 

100 

39 

35 

-180 


GGA+U 5.202 

4.999 

0.753 

0.003 

99 

100 

98 

69 

-140 

Co 2 FeGa 

GGA 5.149 

5.016 

0.885 

0.189 

65 

98 

181 

56 

-62 


GGA+U 5.259 

5.043 

0.772 

0.515 

20 

100 

89 

67 

-151 

Co 2 Fe Alo. 5 Sio. 5 

GGA 5.523 

5.376 

0.755 

0.399 

31 

92 

124 

74 

-119 


GGA+U 5.700 

5.498 

0.667 

0.001 

100 

100 

139 

87 

-125 

Co 2 FeSi 

GGA 5.688 

5.541 

0.714 

2.476 

-55 

86 

189 

24 

-25 


GGA+U 6.196 

5.998 

0.587 

0.008 

98 

100 

73 

54 

-148 

Co 2 FeGe 

GGA 5.854 

5.693 

0.785 

2.288 

-49 

89 

119 

-29 

-49 


GGA+U 6.222 

5.997 

0.624 

0.003 

99 

100 

14 

40 

-570 


P D is positive or negative can be tested by spin-polarized 
angle-resolved photoenrission experiments which, unfor¬ 
tunately, have not been reported on any Heusler com¬ 
pound studied here. 

Table III also indicates that including the on-site 
Coulomb U changes the calculated AHC and SHC sub¬ 
stantially. In particular, the a£ y gets reduced signifi¬ 
cantly for all the studied compounds except Co 2 FeAl (Ta¬ 
ble III). For example, the theoretical <j£ y for Co 2 FeGe is 
119 S/cm from the GGA calculation but is reduced to 
14 S/cm when the on-site Coulomb U is included. This 
suggests that by comparing the calculated a^ y with the 
measured one, one could assess whether or not includ¬ 
ing on-site Coulomb U is needed to properly describe 
the electronic properties of a Co 2 Fe-based Heusler com¬ 
pound. The measured a^ y of Co 2 FeSi 6 is ~160 S/cm, 
being in good agreement with the GGA result (Table 
II). However, it is two times larger than the result of the 
GGA+/7 calculation (about 70 S/cm). This indicates 
that Co 2 Fe-based Heusler compounds are not strongly 
correlated systems and there may be no need to include 
the on-site Coulomb U for these compounds. 


IV. CONCLUSIONS 

We have carried out a systematic ab initio study of 
the anomalous Hall effect and current spin polarization 
as well as the magnetic properties of the Co 2 FeX (X 
= Al, Ga, In, Si, Ge, Sn) Heusler compounds in the 
cubic L2i structure by using the highly accurate all- 
electron FLAPW method. First, we find that the spin- 
polarization of the longitudinal current (P L ) in Co 2 FeX 


(X = Al, Ga and Alo.sSio.s) is ~100 % even though the 
static spin polarization ( P D ) is not. Furthermore, the 
other compounds also have a high current spin polariza¬ 
tion with P L > 85 %. This indicates that all the Co 2 FeX 
compounds are promising for spintronic devices. Inter¬ 
estingly, P D is negative in Co 2 FeX (X = Si, Ge and Sn), 
differing in sign from the P L as well as from that from the 
transport experiments. Second, the calculated AHCs are 
moderate, being within 200 S/cm, and agree well with 
the available experiments on highly L2i ordered Co 2 FeSi 
specimen although they differ significantly from the re¬ 
ported experiments on other compounds where the B2 
antisite disorders were present. Surprisingly, the AHC 
in Co 2 FeSi decreases and then changes from the nega¬ 
tive to positive when Si is replaced by Ge and finally by 
Sn. Third, the calculated total magnetic moments are 
in good agreement with the experiments in all the stud¬ 
ied compounds except Co 2 FeSi where a difference of 0.3 
^b/ f.u. exists. We have also performed the GGA+U 
calculations in order to examine the effects of the on¬ 
site Couloumb repulsion. We find that including the U 
changes the calculated total magnetic moment, spin po¬ 
larization and AHC significantly. In most cases, unfortu¬ 
nately, this results in a worse agreement with the avail¬ 
able experimental results. These interesting findings are 
analyzed in terms of the underlying band structures. 
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